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Abstract 

Data from HERA, LEP and the Tevatron, as well as from low energy experi- 
ments are used to constrain the Yukawa couplings for scalar and vector leptoquarks 
in the Buchmuller-Riickl-Wyler effective model. In the limit of very high lepto- 
quark masses constraints on the coupling to the mass ratio \/M are derived using 
the contact-interaction approximation. For finite masses the coupling limits are 
studied as a function of the leptoquark mass. Some leptoquark models are found 
to describe the existing experimental data much better than the Standard Model. 
Increase in the global probability observed for models including Si or V Q leptoquark 
production/exchange corresponds to more than a 3a effect. Assuming that a real 
leptoquark signal is observed, calculated is an allowed region in the A — M plane. 
The leptoquark signal is mostly resulting from the new data on the atomic parity 
violation in cesium, but is also supported by recent LEP2 measurements, unitarity 
violation in the CKM matrix and HERA high-Q 2 results. 



1 Introduction 

In 1997 the HI |lj and ZEUS || experiments at HERA reported an excess of events in 
positron-proton Neutral Current Deep Inelastic Scattering (NC DIS) at very high mo- 
mentum transfer scales Q 2 , as compared with the predictions of the Standard Model. As 



a possible sign of some "new physics" this results provoked many theoretical specula- 
tions. Clustering of HI events at a positron-jet invariant masses of about 200 GeV was 
considered to indicate a possible resonant production of leptoquark states. The agree- 
ment with the Standard Model prediction improved after both experiments doubled their 
positron-proton data samples, but some discrepancy is still there and calls for a better 
understanding. 

In a paper presented last year ||,[] data from HERA, LEP and the Tevatron, as well 
as from low energy experiments were used to constrain the mass scale of the possible 
new electron-quark contact interactions. A contact interaction model was used as the 
most general framework which can describe possible low energy effects coming from "new 
physics" at much higher energy scales. This includes the possible existence of second- 
generation heavy weak bosons, leptoquarks, as well as electron and quark compositeness 
|5], [J. In addition to the general models, in which all new contact interaction couplings 
can vary independently, the global analysis considered also a set of one-parameter models 
which assumed fixed relations between couplings. However, only parity conserving models 
were selected, as suggested by ZEUS , to avoid strong limits coming from atomic parity 
violation (APV) measurements ||. No significant improvement in the description of the 
data has been obtained for any of these models. 

Theoretical uncertainties in the parity violation measurements in cesium atoms have 
been recently significantly reduced. As a result, measured value of the cesium weak charge 
is now more than 2a away from the Standard Model predictions ||. This discrepancy 
could be due to new parity- violating electron-quark interactions. Considered in this paper 
are effects induced by the possible existence of the first-generation leptoquarks. Predic- 



tions based on the Buchmuller-Riickl-Wyler (BRW) effective model [10| are compared 
with the existing experimental data. In the limit of very high masses, the exchange of 
leptoquarks can be described using the contact interaction approach [n]]. Limits on the 
ratio of the coupling and the mass are derived. For finite leptoquark masses limits on 
leptoquark Yukawa coupling A are studied as a function of the leptoquark mass. 

The aim of the present analysis is to combine the APV measurements with other data 
to constrain leptoquark coupling and mass, and to look for a possible leptoquark signal 
in the combined data. The BRW model used in this analysis is described in section |2|. 
In section R| the relevant data from HERA, LEP, the Tevatron and other experiments are 



1 For recent update of presented results see 
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briefly described. Methods used to compare data with leptoquark model predictions and 
to derive coupling limits are summarised in section |]. The analysis results for different 
leptoquark types, including extracted coupling-mass limits and discussion of the possible 
leptoquark signal are presented in section |[ 

The analysis presented here is based on the approach used in the global analysis of 
eeqq contact interactions || M, which in turn followed [|1^, |13|| . When finalising this 
analysis another work discussing leptoquark exchange as a possible explanation for the 
APV result was released ||14|| . However, the analysis presented there is limited to the 
contact interaction approximation. 



2 Leptoquark models 

Striking symmetry between quarks and leptons in the Standard Model strongly suggests 
that, if there exist a more fundamental theory it should also introduce a more funda- 
mental relation between them. Such lepton-quark "unification" is achieved for example 
in different theories of grand unification (|Tj| and in compositeness models. Whenever 
quarks and leptons are allowed to couple directly to each other, a quark-lepton bound 
state can also exist. Such particles, called leptoquarks, carry both colour and fractional 
electric charge and a lepton number. Also supersymmetric theories with broken R-parity 
predict squarks (leptoquark type objects) coupling to quark-lepton pairs. 

In this paper a general classification of leptoquark states proposed by Buchmiiller, 



Riickl and Wyler [1C] will be used. The Buchmiiller-Ruckl-Wyler (BRW) model is based 
on the assumption that new interactions should respect the SU(3)c x SU{2) L x U{l)y 
symmetry of the Standard Model. In addition leptoquark couplings are assumed to be 
family diagonal (to avoid FCNC processes) and to conserve lepton and baryon numbers 
(to avoid rapid proton decay). Taking into account very strong bounds from rare decays 
it is also assumed that leptoquarks couple either to left- or to right-handed leptons. With 
all these assumptions there are 14 possible states (isospin singlets or multiplets) of scalar 
and vector leptoquarks. Table [I| lists these states according to the so-called Aachen 



notation |16|. An S(V) denotes a scalar (vector) leptoquark and the subscript denotes the 
weak isospin. When the leptoquark can couple to both right- and left-handed leptons, an 
additional superscript indicates the lepton chirality. A tilde is introduced to differentiate 
between leptoquarks with different hypercharge. Listed in Table [i] are the leptoquark 
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Table 1: A general classification of leptoquark states in the Buchmuller-Rtickl-Wyler 
model. Listed are the leptoquark fermion number, F, electric charge, Q (in units of 
elementary charge), the branching ratio to electron-quark (or electron-antiquark), f3 and 
the flavours of the coupled lepton-quark pairs. Also shown are possible squark assignments 
to the leptoquark states in the minimal supersymmetric theories with broken R-parity. 
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fermion number F, electric charge Q, and the branching ratio to an electron-quark pair 
(or electron-antiquark pair), (3. The leptoquark branching fractions are predicted by the 
BRW model and are either 1, | or 0. For a given electron-quark branching ratio (3, the 
branching ratio to the neutrino-quark is by definition (l — (3). Also included in Table [l] are 
the flavours and chiralities of the lepton-quark pairs coupling to a given leptoquark type. 
In three cases the squark flavours (in supersymmetric theories with broken R-parity) 
with corresponding couplings are also indicated. Present analysis takes into account 
only leptoquarks which couple to the first-generation leptons (e, v e ) and first-generation 
quarks (u, d), as most of the existing experimental data constrain this type of couplings. 
Second- and third-generation leptoquarks as well as generation-mixing leptoquarks will 
not be considered in this paper. It is also assumed that one of the leptoquark types 
gives the dominant contribution, as compared with other leptoquark states and that the 
interference between different leptoquark states can be neglected. Using this simplifying 
assumption, different leptoquark types can be considered separately. Finally, it is assumed 
that different leptoquark states within isospin doublets and triplets have the same mass. 

The ep collider HERA is the unique place to search for the first-generation leptoquarks, 
as single leptoquarks can directly be produced in electron-quark interactions. The influ- 
ence of the leptoquark production or exchange on the ep NC DIS cross-section can be 
described as an additional term in the tree level eq —>■ eq scattering amplitude:^ 



l] [S > hU) + sin 2 9 W cos* 9 W t-Ml + ^ lS ' U) 



where s, t and u are the Mandelstam variables describing the electron-quark scattering 
subprocess, e q is the electric charge of the quark in units of the elementary charge, the 
subscripts i and j label the chiralities of the initial lepton and quark, respectively (i, j = 
L, R), and gf and are electroweak couplings of the electron and the quark. In the limit 
Mlq 3> \/s the leptoquark contribution to the scattering amplitude given by ^(s, u) 
does not depend on the process kinematics and can be written as 

2 



M LQ , 



Vi! = <l\£r-\ , (2) 



where M L q is the leptoquark mass, Xlq the leptoquark-electron-quark Yukawa coupling 

and the coefficients are given in Table |2] |TT| . The effect of heavy leptoquark production 

2 Amplitude given for electron-quark scattering describes also scattering of positrons and anti-quarks 
taken with opposite chiralities. 
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Table 2: Coefficients a^J denning the effective contact interaction couplings = at] -j^f- 

LQ 

for different models of scalar (upper part of the table) and vector (lower part) leptoquarks. 
Empty places in the table correspond to = 0. 

or exchange is equivalent to a vector type eeqq contact interaction. It is interesting to 
notice that 5 scalar leptoquark types (S^ , 5*0, Sh 2 , S^ 2 and S1/2) correspond to the same 
contact interaction coupling structures (but opposite coupling signs) as 5 vector models 
(V , V a R , V 1/2 , Vy 2 and Vy 2 respectively). 

For leptoquark masses comparable with the available ep center-of-mass energy it- 
channel leptoquark exchange process and the s-channel leptoquark production have to 
be considered separately. Corresponding diagrams for F=0 and F=2 leptoquarks are 
shown in Figure [I]. The leptoquark contribution to the scattering amplitude can be now 
described by the following formulae: 

• for w-channel leptoquark exchange ( F=0 leptoquark in e~q or e + q scattering, or 
|F|=2 leptoquark in e + q or e~q scattering) 

n eq \ 2 
„f?f s u ) = a v Al ® 

Ml Q -u' 

• for s-channel leptoquark production (F=0 leptoquark in e + q or e~q scattering, or 
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Standard Model F=0 Leptoquark F=2 Leptoquark 



HERA 




Figure 1: Diagrams describing leading order Standard Model processes and leptoquark 
contributions coming from F=0 and F=2 leptoquarks, for NC DIS at HERA, quark-pair 
production cross-section at LEP and Drell-Yan process at the Tevatron, as indicated in 
the plot. 
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|F|=2 leptoquark in e q or e + q scattering) 



a ij A LQ 



M 2 LQ 



is 



M L q 



where Y L q is the total leptoquark width. The partial decay width for every decay 
channel is given by the formula: 



LQ 



8tt(J + 2) 



where J is the leptoquark spin. 



For processes such as e + e — > hadrons a corresponding formula can be written for the 
e + e~ — > qq tree level amplitude: 



4:7C(y ern Cq 



+ 



9-9] 



sin 6w cos 2 9 



(3) 



where the subscripts i and j label the chiralities of the initial lepton and final quark 
respectively and 



n eq \ 2 
a ij A LQ 

M 2 LQ -t 



for F = 



n eq \ 2 

a v JfQ_ for |p. 



u 



Same formulae apply also to qq — > l + l~ amplitude, with % and j labelling the chiralities 
of the initial quark and final lepton respectively. 

Leptoquark states with f3 = | (coupling to both electron-quark and neutrino-quark 
pairs) contribute also to the charged current DIS at HERA eq —> vq' . For M L q ^> yfs 
the effective charged current contact interaction coupling is given by 

v 2 



V 



cc 



V 



euvd 



A 



LQ 



M 



(4) 



3 Experimental Data 
3.1 High-Q 2 DIS at HERA 



Used in this analysis are the 1994-97 data on high-Q 2 e + p NC DIS from both HI | 17|| 
and ZEUS p^j , as well as the recent results from e~p NC DIS scattering [19|, p0| . The 



S 



analysis takes into account expected and measured numbers of events in bins of Q 2 . For 
simplicity let us consider a single Q 2 bin ranging from Q 2 min to Q 2 max - Assume that n SM 
events are expected from the Standard Model. 

The leading-order doubly-differential cross-section for positron-proton NC DIS (e + p — > 
e + X) can be written as 

j2 LO i 

= T ^E^Q 2 ){l^| 2 + |M^| 2 + (l-y) 2 [|M21| 2 + |M^| 2 ]} + 

q(x,Q 2 ) [\MZ\ 2 + \M% R \ 2 + (1 - yf [|M^| 2 + |MfJ 2 ]} , 

where x is the Bjorken variable, describing the fraction of the proton momentum carried 
by the struck quark (antiquark), q(x,Q 2 ) and q(x,Q 2 ) are the quark and antiquark mo- 
mentum distribution functions in the proton and are the scattering amplitudes of 
equation (H), which can include contributions from leptoquark production or exchange 
processes. 

The cross-section integrated over the x and Q 2 range of an experimental Q 2 bin is 



Qmax 1 



d 2 a L °(\ LQ ,M LQ ) 



^ U (X LQ ,M LQ ) = I dQ 2 I dx d J Q2 ^ (5) 



Q 



2 

in i n 



where y ma x is an upper limit on the reconstructed Bjorken variable y, y = — , imposed in 
the analysis. The number of events expected from the Standard Model with leptoquark 
contributions can now be calculated as: 

n(X LQ ,M LQ ) = n SM -( aLO{XL ^ MLQ) ) , (6) 



where is the Standard Model cross-section calculated with formula (|j) (setting A = 0). 
Leading-order expectations of the leptoquark models are used to rescale the Standard 
Model prediction nsM coming from detailed experiment simulation. This accounts for 
different experimental effects, and (to some extent) for higher order QCD and electroweak 
corrections.^ NLO QCD corrections to the resonant leptoquark production are introduced 



as an additional correction factor, based on 121 



For models with leptoquarks coupling to both electron-quark and neutrino-quark pairs 
(S' L , Si, V L and V t ), HERA data on e+p and e'p CC DIS |T7], 0, |2§ are also included 
in the fit. 



3 Correctly taken into account are only those corrections which are the same or similar for the Standard 
Model and for the cross-section including leptoquark contributions. 
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In the limit of heavy leptoquark masses [Mlq 3> \[$) the Q 2 distribution of NC and 
CC DIS events is most sensitive to the leptoquark couplings. For masses below s/s ~300 
GeV, where direct leptoquark production becomes possible at HERA, better limits are 
obtained from studying the electron-jet invariant mass distribution. However, to describe 
correctly the narrow leptoquark resonance production and reconstruction, sizable QED 
and QCD corrections as well as complicated detector effects have to be taken into account. 
As these corrections could not be included in the analysis, the Q 2 distribution was used 
to constrain leptoquark couplings in the whole mass range. Comparison between limits 
calculated from the Q 2 distribution of the ZEUS e + p NC DIS data |18| and the published 
ZEUS limits for F=0 leptoquarks |23[f] is presented in Figure j|. Taking into account 
that ZEUS analysis includes mass dependent selection cuts and that it was optimised for 
leptoquark search, the difference between the two approaches is surprisingly small. Direct 
ZEUS limits are at most 40% lower (depending on the model and the mass range) than 
the one obtained from the Q 2 distribution. 



3.2 Measurements from LEP 

Many measurements at LEP are sensitive to different kinds of "new physics". The lep- 
toquark exchange contribution can be directly tested in the measurement of the total 
hadronic cross-section above the Z° pole.^ The leading order formula for the total quark 
pair production cross-section, a(e + e~ — > qq), at an electron-positron center-of-mass en- 
ergy squared, s, is 



a LO (s) ?>S 



1287T (/ 



J dcos8 



Mrj 2 +|M^|^(l + cos^+(|Ml^ + |M-^(l-cos^j , (7) 

where M^ e are the scattering amplitudes described by equation (|3|), including contribu- 
tions from leptoquark exchange and 9 is the quark production angle in the e + e~ center-of- 
mass system. For comparison with measured experimental values, the expected Standard 
Model cross-section a SM (s) quoted by experiments are rescaled using the ratio of the 



4 Similar limits on the leptoquark couplings and masses have also been presented by the HI Collabo- 
ration |||. 

5 For the leptoquark masses and couplings considered here the effects of the possible leptoquark ex- 
change at \fs = Mz are completely negligible in comparison with the resonant Z° production. 
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Figure 2: Comparison between limits calculated from the Q 2 distribution of the ZEUS 



e + p NC DIS data (this analysis) and the published ZEUS limits [23] for selected F=0 
leptoquarks, as indicated in the plot. 
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leading order cross-sections with and without leptoquark contribution: 

a(s,X LQ ,M LQ ) = a^(s)-( aL ° {sA JS ,MLQ) ) . W 

\ a SM ) 

where a^li is the leading-order Standard Model cross-section (A = 0), calculated with 
equation (0). This takes into account possible experimental effects and higher order QCD 
and electroweak corrections. Included in the analysis are data on <Jhad from Aleph, Delphi, 
L3 and Opal experiments for center-of-mass energies up to 202 GeV f2~5| , |26| , |27| , [28f . All 
measurements are in good agreement with the Standard Model predictions. However, 
cross-section values obtained for 1/^=192-202 GeV are on average about 2.5% above the 
predictions. The combined significance of this deviation is only about 2.3a but it has 
an important influence on the global analysis. 

In the global analysis of electron-quark contact interactions [§, the strongest con- 
straints on the contact interaction couplings resulted from the LEP data on heavy quark 
production, R q (q = b,c), and on forward-backward asymmetries, A q FB . However, this 
is only the case for models assuming family universality. For the first-generation lepto- 
quarks, the constraints resulting from LEP measurements based on heavy flavour tagging 
are much weaker than those resulting from hadronic cross-section measurements. Never- 
theless, possible deviations in the uu and dd quark pair production cross-sections (resulting 
in the deviation of the total hadronic cross-section) can be also constrained using results 
on R c and Rb- Results on A q FB are included in the presented analysis for consistency with 
the previous study 0. 



3.3 Drell-Yan lepton pair production at the Tevatron 

Used in this analysis are data on Drell-Yan electron pair production (pp — > e + e~ X) from 



the CDF |K| an d D0 [31| experiments. The leading order cross-section for lepton pair 
production in pp collisions is 



d 2 a LO Ml 



dM n dY 192tts , ; 

/dcostf [(\Mll\ 2 + \M e R %\ 2 ) (1 + costf) 2 + (|Mf R | 2 + |M- | 2 ) (1 - cos#) s 

where Mu is the invariant lepton pair mass, Y is the rapidity of the lepton pair, 9 is the 
lepton production angle in their center-of-mass system, and %\ and x 2 are the fractions 
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of the proton and antiproton momenta carried by the annihilating qq. When integrating 
over 9, the angular detector coverage is taken into account. The scattering amplitudes 
M^ e and the parton density functions are calculated at a mass scale 

/i 2 = = S = XiX 2 S , 

where s is the total proton-antiproton center-of-mass energy squared. 

The cross-section corresponding to the Mu range from M min to M max is calculated as 

where Y max is the upper limit on the rapidity of the produced lepton pair: 

Y = ln^ 

1 max A/T 

The cross-section calculated with equation @ is used to calculate the number of events 
expected from the Standard Model with leptoquark contribution using formula (0). 



3.4 Direct limits from the Tevatron 

The D0 and CDF experiments at the Tevatron presented limits on the first-generation 
scalar leptoquark masses from the search for leptoquark pair production in hard inter- 
actions (pp — > LQ LQ X). Both experiments see no leptoquark candidate events, with 
leptoquarks decaying into an electron and a jet, above a reconstructed leptoquark mass 
of 200 GeV |32j, The result of the NLO cross-section calculations^ |34| can be param- 
eterised in this mass region as 



The expected number of leptoquark events reconstructed in eejj channel is 



T^exp 



(M LQ ) = eC ■ a s LQ (M LQ ) ^/5 2 Q , 

LQ 



where the sum is over leptoquark states within the considered multiplet and the combined 
effective luminosity (i.e. luminosity corrected for selection efficiency) for two experiments 



3 assuming mass scale /x = 2Mlq 
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is eC « 78 pb" 1 . For leptoquark states with (3 = 0.5, the results of D0 search in eujj 
channel are also included in the analysis. Because of the assumed mass degeneration the 
mass limits for scalar leptoquark multiplets can be significantly higher than for single 
leptoquarks. For Sy 2 isospin doublet (J2 Plq = ^) the combined limit is Mlq > 263 GeV, 



as compared with the published limit of 242 GeV for single leptoquark production [35 



For vector leptoquarks, pair production cross-section at the Tevatron strongly depends 
on the unknown (not constrained in the BRW model) anomalous leptoquark couplings. 
For the presented analysis the values giving the smallest vector leptoquark pair production 



cross-section were assumed |36]. The LO vector leptoquark production cross-section has 
been parametrised as: 

al Q (M LQ ) « 268.1 pb • exp . 

Only the D0 experiment has presented limits on the vector leptoquark masses in the 
minimum cross-section model The limit Mlq > 245 GeV (for (3=1) corresponds to 
the effective luminosity of eC ~ 58 pb -1 . 



3.5 Data from low energy experiments 

The low energy data are included in the present analysis in exactly the same way as in 
the contact interaction analysis j3]. For all leptoquark models the following constraints 
from low energy experiments are considered: 

• Atomic Parity Violation (APV) 

The Standard Model predicts parity non-conservation in atoms caused (in lowest 
order) by the Z° exchange between electrons and quarks in the nucleus. Experi- 
mental results on parity violation in atoms are given in terms of the weak charge 
Qw of the nuclei. Standard Model prediction for Q w are based on the very precise 
measurement of the sin 2 Qw at LEP1 and SLD. A new determination of Qw for 
Cesium atoms was recently reported M. The experimental result differs from the 
Standard Model prediction by: 

AQ% = Q™ S -Q™ = 1.13 ±0.46 

As already mentioned in the Introduction, this 2.5<r discrepancy between the mea- 
surement and Standard Model predictions induces significant evidence for some lep- 
toquark models. Also other "new physics" processes, as for example Z°' exchange, 
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were proposed as a possible explanation for the APV measurement. One has to 
take into account that these new processes can also affect precision measurements 
at LEP1 and the determination of sin 2 Qw, making the analysis much more difficult. 
However, for the leptoquark masses and couplings considered here the effects of the 
possible leptoquark exchange at yfs = Mz can be safely neglected. 

The leptoquark contributions to Qw is: 

97 _|_ AT 

&QW = ^^(V e L U L + V e L U n-VRL-VR R ) 

where rf^ are the effective couplings given by (§). 
Electron-nucleus scattering 

The limits on possible leptoquark contributions to electron-nucleus scattering at low 
energies can be extracted from the polarisation asymmetry measurement 

_ da R - da L 
d<jR + daL 

where daum denotes the differential cross-section of left- (right-) handed electron 
scattering. Polarisation asymmetry directly measures the parity violation resulting 
from the interference between Z° and 7 scattering amplitudes. For isoscalar targets, 
taking into account valence quark contributions only, the polarisation asymmetry 
for elastic electron scattering is 



.4 



3V2G F Q 



2 



2 {gl + g u R ) - (pi + g d R ) 



20na em 

where Q 2 is the four-momentum transfer and the effective electroweak coupling of 
the quark is modified by the leptoquark contribution 



9l 



eq 

q VLi / -1 «-| \ 

eff 9i ~ 2V2G F ' 



The data used in this analysis come from the SLAC eD experiment pg[ , the Bates 
eC experiment |39| and the Mainz experiment on eBe scattering [^Oj] . 



For leptoquarks contributing to charged current processes, additional constraints come 
from: 
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Lepton-hadron universality of weak Charged Currents 

New charged current interactions would affect the measurement of V u d element of 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix, leading to an effective violation 



of unitarity |41| , |42fl . The new experimental constraint is [[43 



\V ud \ 2 + \V US \ 2 + \V ub \ 2 = 0.9959 ±0.0019 



whereas the expected leptoquark contribution is 



v ud = v™- i " 



with 7] is given by equation ([|). 
Electron-muon universality 

In the similar way new charged current interactions would also lead to effective 
violation of e-fi universality in charged pion decay j|l|. The current experimental 
value of R = T(it~ — > ev)/T{it~ — > fiu) is fl44[ 



R 

-p^ = 0.9966 ±0.030 , 

Rsm 

whereas the expected contribution from leptoquark exchange is 



R 



Rsm • ( 1 ' 



2 v / 2G i 



It is important to notice, that data in the charged current sector also indicate a possible 
deviation from the Standard Model predictions: a slight violation of the unitarity of the 
CKM matrix and of the e-fi universality. The combined significance of these two results 
is about 2. 4a and has a considerable influence on the presented analysis. 



4 Analysis method 

The analysis method is similar to the one used in the recently published analysis 0. 
For every leptoquark coupling and mass value the probability function describing the 
agreement between the model and the data is calculated: 

V(X LQ ,M LQ ) ~ WP^Xlq.Mlq). (11) 



16 



The product runs over all experimental data i. The logarithm of the probability function 
InV is the so called log-likelihood function, which is often used in similar analysis: 



\nV(X LQ ,M LQ ) = Y,^ P ^ X lq,M lq ). 

i 

The data used in this analysis can be divided into two classes. 

1. For experiments in which a result is presented as a single number with an error 
which is considered to reflect a Gaussian probability distribution, the probability 
function can be written as 

\2\ 



Pi(^LQ, M L q) ~ exp -- v v ^ '-] , (12) 



\{F{\ LQ ,M m )-AAf 
A / 

where AA is the difference between the measured value and the Standard Model 
prediction, a a is the measurement error and F(X L q, M L q) is the expected leptoquark 
contribution to the measured value. This approach is used for all low energy data 
as well as for the LEP hadronic cross-section measurements. 

2. On the other hand, if the experimentally measured quantity is the number of events 
of a particular kind (e.g. HERA high-Q 2 data or Drell-Yan lepton pairs and direct 
search results from the Tevatron), and especially when this number is small, the 
probability is better described by the Poisson distribution 

p /x w \ n(A LQ , M LQ ) N ■ exp(-n(\ LQ , M LQ )) 
Fi{\ LQ ,M LQ ) ~ — , {IS) 

where N and ti(Xlq, Mlq) are the measured and expected number of events in a 
given experiment, respectively, and u^Xlq, Mlq) takes into account a possible lep- 
toquark contribution. This approach has been used for the HERA and the Tevatron 
data. 



For low energy data the total measurement error can be used in (|T2"| ) taking into 
account both statistical and systematic errors. For collider data, formula (|T^) or (IB) is 



used to take into account the statistical error of the measurement only. The systematic 
errors are assumed to be correlated to 100% within a given data set (e.g. e + p NC DIS data 
from ZEUS ) they This approach, as well as the migration corrections used for HERA 
and Tevatron Drell-Yan results are discussed in detail in B. 



17 



The probability function V(Xlq, Mlq) summarises our current experimental knowl- 
edge about possible leptoquark couplings and masses. As V is not a probability distribu- 
tion, it does not satisfy any normalisation condition. Instead it is convenient to rescale 
the probability function in such a way that for the Standard Model it has the value of 1: 

V(X LQ = 0,M LQ ) = 1. (14) 
or \nV(X LQ = 0,M LQ ) = 0. 

Using the probability function V(\lq, Mlq) two types of limits in (Xlq, Mlq) space 
are calculated: 

• Rejected are all models (parameter values) which result in 

V(X LQ ,M LQ ) < 0.05 
or lnV(X LQ ,M LQ ) < -3.0 

This is taken as the definition of the 95% confidence level (CL) exclusion limit. 
Exclusion limits presented in this paper are lower limits in case of leptoquark mass 
Mlq and upper limits in case of Xlq or Xlq /Mlq. 

• Some leptoquark models turn out to describe the data much better than the Stan- 
dard Model: 

Vmax = max V(X L q,M L q) > 1. 

In that case the 95% CL signal limit corresponding to the uncertainty on the 
"best" values of Xlq and Mlq is defined by the condition 

v(x L q,m L q) > om-v max 

or \iiV(Xlq,Mlq) < \nV max - 3.0 

In the previous analysis Q no significant deviations from the Standard Model were 
observed. In such a case both definitions give similar results and there is no need to 
distinguish between exclusion and signal limits. 
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Table 3: Coupling to mass ratio, (tt 2 - J , resulting in the best description of the ex- 

\ 1V1 LQ / max 

perimental data in the contact interaction approximation, and the corresponding model 
probability V max and the log-likelihood lnV max , for different leptoquark models, as in- 
dicated in the table. The errors attributed to non-zero values correspond to the 
decrease of \nV by |. Also given are 95% CL signal (for models with V ma x > 20) and 
(upper) exclusion limits on jjj^, and (lower) exclusion limits on leptoquark masses Mlq. 

5 Results 

In the limit of very high leptoquark masses (contact interaction approximation) the prob- 
ability function depends only on the Xlq/M L q ratio. Using the global model probability 
V(X LQ , M L q), as defined by equation (||), the value {X LQ /M LQ ) max giving the maximum 
probability is determined for each model. The results are presented in Table |3]. The 
attributed errors, quoted for models which give better description of the data than the 
Standard Model (i.e. (Xlq/Mlq) ^ > 0) correspond to the decrease in lnP(Ax,Q, M L q) 
by |. The probability functions V(Xlq, Mlq) for different leptoquark models are shown 
in Figure |[ 

For 8 out of 14 leptoquark models, the Standard Model gives the best description 
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Figure 3: Probability function V(\lq, Mlq) (left hand scale) and the log-likelihood 
function 1i[iV(\lq, M^q) (right hand scale), in the limit of very high leptoquark masses, 
for different leptoquark models as indicated on the plot. 
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of the considered experimental data ((Xlq/ '^LQ) max = 0). 95% CL exclusion limits for 
Xlq/M L q range for these models from 0.12 TeV -1 (for V Q L model) to 0.29 TeV -1 (for Sy 2 
model). The other 6 models are able to describe the data better than the Standard Model. 
In all cases the "best" coupling to mass ratio turns out to be of the order of 0.3 TeV -1 . 

The best description of the data is given by the Si model for (Xlq/Mlq)^^,, = 
0.28 ± 0.04 TeV -1 resulting in the maximum probability V max =367 (lnP maa; =5.9). For 
the Gaussian probability function this would correspond to about 3.4a deviation from 
the Standard Model. The effect is mainly due to the APV result: the contribution of the 
APV measurement to the maximum probability is P=20 (lnP=3.0), corresponding to a 
2.4a deviation from the Standard Model. The result is also supported by the low energy 
charged current data (unitarity of the CKM matrix and e-/i universality; lnP=2.4, 2.2cr 
effect) and LEP2 hadronic cross-section measurements (lnP=0.5, 1.0a effect). Contribu- 
tions of different data sets to the Si model probability function are presented in Figure |J 
The fitted value of (Xlq/ '^LQ) max results in almost the best description of both APV and 
low energy charged current data, whereas LEP2 hadronic cross-section measurements sug- 
gest even higher values of Xlq/Mlq ~ 0.7 TeV -1 . The 95% CL signal limit, corresponding 
to model the probabilities V > 0.05 • V max is 0.15 < X LQ /M LQ <0.36 TeV -1 . 

The V model also gives a very good description of the data, resulting in V max =Yl1 
(lnP=4.8 corresponding to about 3.1<r). In this case the APV result (lnP=3.1, 2.5a) is 
strongly supported by LEP2 data (lnP=1.3, 1.6a). The Sy 2 and Vy 2 models describe 
the APV measurement as well but they do not improve the description of other data. For 
and Vi/2 models, the coupling values required to explain APV data are disfavoured 
by other experiments (mainly by LEP2 hadronic cross-section measurements) resulting 
in even smaller V ma x values. Signal limits for 4 models which result in V max > 20 are 
included in Table |3|. For models with V max > 1 (models describing the APV data) the 
95% CL exclusion limits on Xlq/M L q range from 0.41 TeV -1 (for Si model) to 0.52 TeV -1 
(for V model). 

All of the results presented above were based on the contact interaction approximation, 
which is valid for the leptoquark masses above about 1 TeV. In the second part of the 
presented analysis lower leptoquark masses were also considered. In that case, leptoquark 
constraints have to be studied in terms of the leptoquark coupling and the leptoquark 
mass as two independent parameters. 

Below 1 TeV, effects of the finite leptoquark mass reduce the virtual leptoquark ex- 
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Figure 4: Contributions of different data sets (as indicated on the plot) to the global 
probability functions V(Xlq, M^q) (left hand scale) and to the log-likelihood function 
\iiV(Xlq, M L q) (right hand scale), for the S 1 model in the limit of very high leptoquark 
masses. 
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change contribution to the expected LEP and Tevatron cross-sections. However, this 
effect is small and the asymptotic limit on Xlq/M L q increases by only about 10% for lep- 
toquark masses Mlq ~ 300 GeV. For masses below 300 GeV, the limits on Xlq become 
much stronger because of the direct searches at HERA and at the Tevatron. Combined 
constraints on the leptoquark coupling and mass are derived from the probability func- 
tion V(\lq, Mlq), as described in section [| The 95% CL exclusion limits in (Xlq, Mlq) 
space, for different models of scalar and vector leptoquarks are presented in Figure The 
95% CL exclusion limits on the leptoquark masses (i.e. largest mass values resulting in 
V < 0.05 for any value of Xlq) are included in Table ||. 

The parameter values resulting in the best description of the experimental data were 
also searched for for finite leptoquark masses, varying \ L q and M L q independently. Only 
for one leptoquark model an improvement has been obtained as compared with the asymp- 
totic solution. For Mlq=276 GeV and Alq=0.095 the maximum probability P ma:r =142 
(hiP=5.0) is obtained for the V model. This corresponds to about 3.1a deviation from 
the Standard Model. The effect is mainly due to the new APV measurement (lnP=3.0, 
2.4a effect), but is also supported by the excess of high-Q 2 NC e + p DIS events at HERA 
(lnP=1.4, 1.7a effect) and the LEP2 hadronic cross-section measurements (lnP=1.2, 
1.5(7 effect). For all HERA, LEP and low energy data the maximum probability turns out 
to be V max =367 (as compared to V max =122 in the contact interaction limit). However, 
the value of V leptoquark mass of M L q=276 GeV is already strongly disfavoured by the 
negative direct search results from the Tevatron (P=0.36, lnP=-1.0). Contributions of 
different data sets to the probability function for the V model with Mlq = 276 GeV 
are presented in Figure |J Very good description of APV and HERA high-Q 2 data is 
obtained for the fitted value of Xlq, whereas LEP2 measurements again suggest higher 
values of Xlq ~ 0.16. The ratio of the predicted e + p cross-section at high Q 2 to the 
Standard Model cross-section is shown in Figure |^ together with the corresponding HI 



1% and ZEUS |L8[] data. The hypothesis of the V leptoquark production can describe 
the excess of events at highest Q 2 not affecting the perfect agreement with the Standard 
Model at Q 2 < 10000 GeV 2 . Also shown in Figure [7| is the predicted deviation of the 
total hadronic cross-section at LEP as a function of ^/s. Best fit of the V model results 
in the cross-section increase at highest a/s by about 1%, which is consistent with available 
data. From the fit of a two-dimensional Gaussian distribution in the close neighbourhood 
of the maximum of the probability function V(Xlq, Mlq), the errors on the V parameter 
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Figure 5: 95% CL exclusion limits in (Xlq, Mlq) space, for different leptoquark models, 
as indicated in the plot. Excluded are coupling and mass values above or to the left of 
the limit curves. 
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Figure 6: Contributions of different data sets (as indicated on the plot) to the global 
probability functions V(Xlq, Mlq) (left hand scale) and to the log-likelihood function 
In V(Xlq, M lq ) (right hand scale), for the V model with M LQ = 276 GeV. 

values were estimated: 

M LQ = 276 ±7 GeV 
\ LQ = 0.095 ±0.015 . 

The local maximum of the probability function at M^q=274 GeV is also observed for the 
Sy 2 model (P=20.3, as compared with V ma x=35.8 obtained in the high leptoquark mass 
limit). This maximum is due to APV and HERA data, but is strongly suppressed by the 
Tevatron direct search results.[] Signal limits in the (Xlq, Mlq) space were studied for 
all leptoquark models which resulted in the description of the experimental data much 
better than the Standard Model (V ma x > 20). Best parameter values and estimated 95% 
CL lower limits on the leptoquark masses are summarised in Table |j. In Figure the 
signal limits at 68% and 95% CL are compared with exclusion limits in the (Xlq, Mlq) 
space. 



7 For the Sy 2 isospin doublet the combined Tevatron 95% CL limit is Mlq > 263 GeV, as compared 



with the published limit of 242 GeV for single leptoquark production (see Section 3.4) 
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Figure 7: Cross-section deviations from the Standard Model resulting from the fit of the 
V model (thick solid line) compared with HERA NC e + p DIS cross-section results (upper 
plot) and LEP2 hadronic cross-section results (lower plot). 
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Figure 8: Signal limits on 68% and 95% CL for different leptoquark models, as indicated 
in the plot. Dashed lines indicate the 95% CL exclusion limits. For V model a star 
indicates the best fit parameters. For other models the best fit is obtained in the contact 
interaction limit Mlq — > oo. 
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Table 4: Coupling XiQ\ max and mass M L Q\ max values resulting in the best description 
of the experimental data and the corresponding model probability V max for different 
leptoquark models, as indicated in the table. Also given are 95% CL lower limits on the 
leptoquark mass (signal limits). Shown in the table are only those models which give 
much better description of the experimental data than the Standard Model (V ma x > 20). 



When the best description is obtained in the very high mass 




is given. 
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6 Summary 



Data from HERA, LEP and the Tevatron, as well as from low energy experiments were 
used to constrain Yukawa couplings for scalar and vector leptoquarks in the Buchmiiller- 
Riickl-Wyler effective model. In the limit of very high leptoquark masses, constraints 
on the coupling to mass ratio were studied using the contact-interaction approximation. 
Some leptoquark models are found to describe the existing experimental data much better 
than the Standard Model. The best description of the data is obtained for the Si model 
with M LQ > 300 GeV and X LQ /M LQ = 0.28 ± 0.04 TeV" 1 and for the V Q model with 
M L q = 276 ± 7 GeV and Xlq = 0.095 ± 0.015. In both cases the increase of the global 
probability corresponds to more than 3a deviation from the Standard Model. The effect 
is mainly due to the new data on atomic parity violation in cesium, but is also supported 
by LEP2 hadronic cross-section results and HERA NC e + p DIS (for the V model) or low 
energy CC data (for the Si model). Other data considered in this analysis are also in 
good agreement with predictions these models. 

If the observed V D signal is real it could become visible in the new HERA e + p data, 
which are now being collected at increased center-of-mass energy.^] 
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8 Since July 1999 HERA scatters 27.5 GeV positrons on 920 GeV protons resulting in the center-of- 
mass energy • v /s=318 GeV. The currently available HERA e + p data were collected 1994-97 with a proton 
beam energy of 820 GeV, corresponding to ^=300 GeV. 
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